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Abstract 

We investigate the critical behavior that d-dimensional systems with short-range 
forces and a n-component order parameter exhibit at Lifshitz points whose wave- 
vector instabihty occurs in an m-dimensional isotropic subspace of M'^. Utihzing 
dimensional regularization and minimal subtraction of poles in d = 4: + Y~^ 
dimensions, we carry out a two-loop renormalization-group (RG) analysis of the 
field-theory models representing the corresponding universality classes. This gives 
the beta function (3u{u) to third order, and the required renormalization factors as 
well as the associated RG exponent functions to second order, in u. The coefficients 
of these series are reduced to m-dependent expressions involving single integrals, 
which for general (not necessarily integer) values of m G (0, 8) can be computed 
numerically, and for special values of m analytically. The e expansions of the critical 
exponents 77^2, rjn-, 1^12, ^li-, the wave- vector exponent /3g, and the correction-to- 
scaling exponent are obtained to order e^. These are used to estimate their values 
for d = 3. The obtained series expansions are shown to encompass both isotropic 
limits m = and m = d. 

Key words: field theory, critical behavior, anisotropic scale invariance, Lifshitz 
point 



1 Introduction 



The modern theory of critical phenomena [1-3] has taught us that the standard 
101^ models with a n-component order-parameter field (f) = {(pi, i = 1, . . . ,n) 
and 0{n) symmetric action have significance which extends far beyond the 
models themselves: They describe the long-distance physics of whole classes 
of microscopically distinct systems near their critical points. In fact, they 
are the simplest continuum models representing the 0{n) universality classes 
of (i-dimensional systems with short-range interactions whose dimensions d 
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exceed the lower critical dimension (= 1 or 2) for the appearance of a 
transition to a phase with long-range order, and are less or equal than the 
upper critical dimension d* = 4 (above which Landau theory yields the correct 
asymptotic critical behavior) . By investigating these models via sophisticated 
field theoretical methods [4] , impressively accurate results have been obtained 
for universal quantities such as critical exponents and universal amplitude 
ratios. 

A well-known crucial feature of these models is their scale (and conformal) in- 
variance at criticality: The order-parameter density behaves under scale trans- 
formations asymptotically as 



in the infrared limit £ — > 0, where = {d — 2 + ri)/2 is the scahng dimension 
of (p. This scale invariance is isotropic inasmuch as all d coordinates of the 
position vector x are rescaled in the same fashion. 

There exists, however, a wealth of phenomena that exhibit scale invariance 
of a more general, anisotropic nature. Roughly speaking, one can identify 
four different categories: (i) static critical behavior in anisotropic equilibrium 
systems such as dipolar-coupled uniaxial ferromagnets [5] or systems with 
Lifshitz points [6,7], (ii) anisotropic critical behavior in stationary states of 
nonequilibrium systems (like those of driven diffusive systems [8] or encoun- 
tered in stochastic surface growth [9]), (iii) dynamic critical phenomena of 
systems near thermal equilibrium [10], and (iv) dynamic critical phenomena 
in nonequilibrium systems [8]. 

In the cases of the first two categories, the coordinates x can be divided into 
two (or more) groups that scale in a different fashion. Writing x — (ajy, ccj.), 
we call these parallel and perpendicular, respectively. Instead of Eq. (1) one 
then has 



where 9, the anisotropy exponent, differs from one. Categories (iii) and (iv) 
involve genuine time-dependent phenomena for which time typically scales 
with a nontrivial power of the length rescaling factor £. For phenomena of 
category (ii), one cannot normally avoid to deal also with the time evolution. 
This is because a fluctuation-dissipation theorem generically does not hold for 
such nonequilibrium systems; their stationary-state distributions are not fixed 
by given Hamiltonians of equilibrium systems and hence have to be determined 
from the long-time limit of their time-dependent distributions in general. 

Category (i) provides very basic examples of systems exhibiting anisotropic 



<t){ix) - ^-^^ (f){x) 




(l>{fx\i,ex^)^r^^ 4>{x\^,x^), 



(2) 
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scale invariance whose advantage is that they can be investigated entirely 
within the framework of equilibrium statistical mechanics. The particular ex- 
ample we shall be concerned with in this paper is the familiar continuum 
model for an m-axial Lifshitz point, defined by the Hamiltonian 

H = ji'x{q + f {A,^f + \ (V.4.f + f 0^ + 1^ I'^r'} ■ (3) 

Here 4>{x) = {4>i{x))^^^ is an n-component order-parameter field where x = 
(a;||,a;_L) G IR"" x R'^-"^. The operators Vy, V_l, and Ay denote the d\\ = m 
and d± — d—m dimensional parallel and perpendicular components of the 
gradient operator V and the associated Laplacian Ay = Vy^, respectively. 
The parameters ctq and uq are assumed to be positive. At zero-loop order 
(Landau theory), the Lifshitz point is located at Po = tq = 0. 

We recall that a Lifshitz point is a critical point where a disordered phase, 
a spatially uniform ordered phase, and a spatially modulated ordered phase 
meet. For further background and extensive lists of references, the reader is 
referred to review articles by Hornreich [6] and Selke [7] and to a number of 
more recent papers [11-18]. 

An attractive feature of the model (3) is that the parameter m can be varied. 
It was studied many years ago [19-22] by means of an e expansion about the 
upper critical dimension 

d*(m)=4 + m/2, m<8. (4) 

The order-e results for the correlation-length exponents vi2 and z/;4, first de- 
rived by Hornreich et al. [19], are generally accepted. Yet long-standing contro- 
versies existed on the terms of the correlation exponents rji2 and ?7;4 and the 
wave- vector exponent (5q. Mukamel [20] gave results for all m with 1 < m < 8. 
These agreed with what Hornreich and Bruce [21] found in the uniaxial case 
111 = 1 via an independent calculation, but were at variance with Sak and 
Crest's [22] for m = 2 and m — Q (who investigated only these special cases). 

More recently, Mergulhao and Carneiro [13,14] presented a reanalysis of the 
problem based on renormalized field theory and dimensional regularization. 
Treating explicitly only the cases m = 2 and m = 6, they recovered Sak 
and Crest's results for 1^12 and ?7;4, but did not compute j3q. They analytically 
continued in dy rather than in and, fixing the latter at = 4 — m/2 while 
taking the former as d\\ — m — ey, with m = 2 or 6, they also derived the 
expansions of the correlation-length exponents V12 and to order ey. 

The purpose of the present paper is to give a full two-loop renormalization 
group (RC) analysis of the model (3) for general, not necessarily integer values 
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of m G (0,8) in d = d*{m) — e dimensions. Q As a result we obtain the e 
expansions of all critical exponents rii2, rji^, Pg, 1^12, and to order e^. In a 
previous paper [18], hereafter referred to as I, we have shown how to overcome 
the severe technical difficulties that had hindered analytical progress in this 
field and prevented a resolution of the above-mentioned controversy for so 
long. Working directly in position space and exploiting the scale invariance of 
the free propagator at the Lifshitz point, we were able to compute the two- 
loop graphs of the two-point vertex function F^^^ and T^^^^^^, its analog with 

an insertion of / d'^a; (V||0)^/2. Together with one-loop results, these suffice for 
determining the exponents rii2, rji^, and Pg to order e^. In order to obtain the 
correlation-length exponents 1^12 and z/^ to this order in e we must compute 
the two-loop graphs of the four-point vertex function F'^'^^ and of F^a''- 

Our results are of importance to recent work on the generalization of confor- 
mal invariance to anisotropic scale invariant systems [25-27]. Some time ago 
Henkel [25] proposed a new set of infinitesimal transformations generalizing 
scale invariance for systems of this kind with an anisotropy exponent 6 = 2/p, 
p = 1, 2, .... He pointed out that the case p = 4, 6* = 1/2, is realized for the 
Lifshitz point of a spherical (n — 00) analog of the ANNNI model [28,29], 
and that the same m-independent value of 6 in Ref. [28] was found to persist 
to first order in e for the Lifshitz point of the model (3). However, as can be 
seen from Eq. (66) below [and Eq. (84) of I], 6 deviates from 1/2 at order e^. 
This shows that the Lie algebra discussed in Ref. [25] cannot strictly apply 
below the upper critical dimension (e > 0) if n is finite, except in the trivial 
Gaussian case uq = 0. 

The remainder of this paper is organized as follows. In the next section we 
recapitulate the scaling form of the free propagator for tq = po = 0. We give 
the explicit form of its scaling function as well as those of similar quantities, 
and discuss their asymptotic behavior for large values of their argument. These 
informations are required in the sequel since the expansion coefficients of our 
results for the renormalization factors and critical exponents can be expressed 
in terms of single integrals involving these functions. 

In Sec. 3 we specify our renormalization procedure and present our two-loop 

^ Another two-loop calculation was recently attempted by de Albuquerque and 
Leite [23,24]. In their evaluation of two-loop graphs — e.g., of the last graph of F^^^ 
shown in Eq. (B.7), — they replaced the integrand of the double momentum inte- 
gral by its value on a line. We fail to see why such a procedure, by which more or 
less arbitrary numbers can be produced, should give meaningful results. Let us also 
emphasize that using such 'approximations' leads to the following problem: Unless 
the corresponding 'approximations' are made for higher-loop graphs involving this 
two-loop graph as a subgraph, pole terms that cannot be absorbed by local coun- 
terterms are expected to remain because they will not be canceled automatically 
through the subtractions provided by counterterms of lower order. 
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results for the renormalization factors. Our e-expansion results for the critical, 
correction-to-scaling, and crossover exponents are described in Sec. 4. Utiliz- 
ing these we determine numerical estimates for the values of these exponents 
in d — 3 dimensions, which we compare with available results from Monte 
Carlo calculations and other sources. Section 5 contains a brief summary and 
concluding remarks. In the Appendixes A-E various calculational details are 
described. 



2 Scaling functions of the free theory and their asymptotic behav- 
ior 



2.1 The free propagator and its scaling function 



Following the strategy utilized in I, we employ in our perturbative renormal- 
ization scheme the free propagator with Tq — po — 0. In position space, it is 
given by 

G(a.) = G(x||,x^)= / ""^l , . (5) 

J q± + cToQw 



Here x\\ = \x\\\ and xj_ = \x^\ are the Euclidean lengths of the parallel and 
perpendicular components of x, and we have introduced the notation 

with and / 1^ (6) 

J J (27r)"^ J J (2nY-"' ^ ' 

qii qx 1\\ K"" ' ' qx Rd-m 



for integrals over momenta q = (^ii,^^) G x M*^ Whenever necessary, 
these integrals are dimensionally regularized. 

Rescaling the momenta as qnCo^^y^^ — > q\\ and q±x± — > q^^ yields the scaling 
form [18] (cf. Ref. [25,29]) 

x^) = <7o-"^/' $(ao-^/' xl'/') (7) 



with the scaling function 



^{v) = ^{v;m,d) = j j 



qii q± 
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where e± is an arbitrary unit d_L-vector while v stands for the dimensionless 
dii-vector 



-1/4 

V — aQ ' x\\x 



-1/2 



(9) 



In I the foUowing representation of ^{v) in terms of generahzed hypergeometric 
functions was obtained: 



^v)=2-^-"'7r-—^v) 



(10) 



r(i-f 



iFsl 1- 



ell m 
2' 2' 2^T' 64 



r(i-f) 



'3 e 3 



4 r(l+f) '^A2 2'2'"^^4'64j ' 



(11) 



with e = 4 + y — (i. Upon expanding the hypergeometric functions in powers 
of and resumming, one arrives at the Taylor expansion 



. _ ~ 1 r(i- l + l) 



(12) 



The result tells us that $ can be written in the form 



2' 2 



1 m 1\ i; 

2^T' 2;'~T 



(13) 



where i\E'i is a particular one of the Fox- Wright functions (or Wright func- 
tions) [30-34], further generalizations of the generalized hypergeometric 
functions pFg whose series representations are given by 

,*,[{(... A,)}, {(6,, = |1 ^r^l^^^ ^' . (14) 



In the sequel, we shall need the asymptotic behavior of $(1;) as t; — > cxd. This 
may be inferred from theorems due to Wright [31,32] about the asymptotic 
expansions of the functions p^^. We discuss this matter in Appendix A, where 
we show that the asymptotic expansion these theorems predict for nonexcep- 
tional values of m and e, 

^, . /t;\-^+2^^ (-1)^ 2r(2-e + 2fc) 
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follows from the integral representation (8) in an equally straightforward man- 
ner as the Taylor expansion (12). Nonexceptional values of m and e are char- 
acterized by the property that none of the poles which the nominator of the 
coefficient 



r(i-f + 



r(i + f + 1) 



(16) 



of the power series (12) has at 

k^ki = e-2l , 1^1,2, 



(17) 



gets canceled by a pole of the denominator. If e = 0, the only values among 
m — 1, 2, . . . , 7 for which such cancellations occur are m — 2 and m — 6. More 
generally, this happens for d = m-|- 1 and d — 111 + 3 where the expansion (15) 

terminates after the first {k = 0) term and vanishes identically, respectively. In 
accordance with Wright's theorems, corrections to these truncated expansions 
are exponentially small. In fact, in these two cases ^{v) reduces to the much 
simpler expressions 



$(v) 



u 



2-m 



8 W2 
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m — 2 V 



' 4 



d — m + 1 , 



(18) 



and 



$(i;) = (47r)- 2 e-^ d = m + 3. 



(19) 



where 7(0, x) is the incomplete gamma function. These equations comprise two 
cases where d becomes the upper critical dimension (4), namely {d, m) — (6, 7) 
and {djTn) = (2,5). In the former, Eq. (18) simplifies to 



(27r)3 



1 + 



m 



6 , d = d* = 7. (20) 



This result as well as Eq. (19) with m — 2 were employed in I, where we also 
derived the leading term {k — 0) of the asymptotic series (15). 

For general values of m and d — d*, the scaling function $(v;m,d) = $(v) 
can be written as 



^{v] m, d* 



_t-t 1 6+m 



l-f 2y^, 2' 4 '64/ 



— y/n 



4 I 4 
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where haiz) and Ia{z) are modified Struve and Bessel functions, respectively 
[35]. 

Tfie second form is in conformity with the one given by Frachebourg and 
Henkel [29] for the case m = 1. These authors encountered this (and similar) 
scaling functions when studying Lifshitz points of order L — 1 of spherical 
models. They also analyzed the large-i; behavior of these functions, verifying 
explicitly the asymptotic forms predicted by Wright's theorems. If we let a = 1 
and set x = v^/A, their scaling function denoted ^'(a, x) corresponds precisely 
to our ^{v), and the asymptotic expansion they found is consistent with ours 
in Eq. (15) and the large-?; form (22) presented below.0 

Finally, let us explicitly give the large-i; forms of ^{v,m,d*) as implied by 
Eq. (15): 



^{v;m,d*) 



\l—m 



TT 



r(^ 



jj+rn m-2 \ 1 24 (m- 6) 
960 (m- 10)(m-6) 

102 



+ 



O V 



-16 



f22) 



In accordance with our previous considerations, all terms or all but the first 
one of this series vanish when m = 2 or m = 6, respectively. 



2.2 Other required scaling functions of the free theory 



Besides $, our results to be given below involve two other scaling functions. 
One is the function S(t;; m, d) of I. This is defined through 

- (V||G* V||G)(x) =xI^+Vo-('"^^)/^H(ao%x:^/^), (23) 

where the asterisk indicates a convolution, i.e., {f*g){x) = Jd'^x'f{x — 
x') g{x'). The explicit form of S(t;; m, d*) may be gleaned from Eq. (A5) of I, 

^ Note that the Hamiltonians of the spherical models considered in Ref. [29] involve 
instead of (A||0)^ a derivative term of the form Yl^i{df (f))'^ , which breaks the 
rotational invariance in the parallel subspace R™. Comparisons with our results for 
the free theory are therefore only possible for L = 2 and m = 1. 
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where it was given in terms of Bessel and hypergeometric functions. This can 
be written more compactly as 



S(v; m, d*) 



9 m 

V 2 



g , m 4+m 




(24) 



The asymptotic expansion of the difference of functions in the square brackets 
of this equations follows from Eqs. (12.2.6) and (9.7.1) of Ref. [35], implying 



'E,{v; m, d*) 



-2-m 



TT 



6+m m — 2 



V 



1 + 



6 — m 4^ 



+ 



6 - m 3 (10 - m) 4 



(26) 



The leading term ~ v""^ was already given in I. Note also that again all terms 
or all but the first one of this series vanish when m = 2 or m = 6, respectively, 
as is borne out by the explicit forms 



1 e-"V4 
(.; 2, 5) = -$(.; 2, 5) = ^ 



(26) 



and 

1 _ e-^V4 

The third scaling function we shall need is defined through the Taylor expan- 
sion 

e(v]m) = y T — ^ ■ (28) 

It can be expressed in terms of generalized hypergeometric functions pFq by 
summing the contributions with even and odd k separately. One finds 



V 



V 



32 r(|+f ) 



4 rri- 



/ 3 3 m V 

1^1, 1; 2^2, 2+ j;^^ 

/l 3 m 

^^n2'2'^+T'64 • 



(29) 
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Details of how this function arises in the computation of the Laurent expansion 
of the four-point graph may be found in Appendix C. In Appendix D we 
show that Q{v,m) behaves as 



Q{v; m) 



r(^) 



- In — + ^ 
lb 



Ce-S 



m — 2 



+ 96 



(m-2)(m-6) 



+ 0{v-'') 



(30) 



in the large-f hmit, where Ce — 0.577216 is Euler's constant and while ip{x) 
denotes the digamma function. 

For the special values m = 2 and m = 6, the function 0(i;;m) reduces to a 
sum of elementary functions and the exponential integral function £"1(0;). As 
can be easily deduced from the series expansion (28), one has 



e{v,2) = -2 



CE + ln- + EJ- 



(31) 



and 



e{v,6)^l-2CE-ln--2EA-j + 
= l + e{v;2)-^v;6,7) , 



+ 32 



(32) 



respectively. 



3 Renormalization 



3. 1 Reparametrizations 



^Prom I we know that the ultraviolet singularities of the A^-point correla- 
tion functions (n.^=i(l)i^{x^)) of the Hamiltonian (3) can be absorbed via 
reparametrizations of the form 

= 0^.^^^ ^ (33) 

To - Toc ^ fl^Z^T , (34) 

ao^ Z^a , (35) 
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UoCTq 



-m/4 TP ^ Z V 
^ m,e — H' ) 



(36) 



and 



(Po - Poc) ctq ^^'^ ^ iiZpp . 



(37) 



Here /i is an arbitrary momentum scale. The critical values of the Lifshitz 
point, (Joc and poc, vanish in our perturbative RG scheme based on dimen- 
sional regularization and the e expansion. The factor Fm^e serves to choose a 
convenient normalization of u. A useful choice is to write the following one- 
loop integral for XZX ^ 



1 



(38) 



This integral is evaluated in Appendix C. The result, given in Eq. (B.13), 
yields 



8+m-2e 

i^m,6 = (47r) 



r(i + t) 






r(f) 


r(2 - 


s)r(f; 


) 



(47r) 



r(f: 



l + (2-CE + ln47r)- + 0(e2) 



(39) 



3.2 Renormalization factors 



With this different choice of normalization of the coupling constant, our results 
of I for Z^, Zfj, and Zp translate into 



_ n+2 j^jm) 3 
^^-^"^12f8-m)T + ^^^^' 



^''^^ - ^ + ^ 96m(m+2) T + ^^"^ ^ ' 



(40) 
(41) 



and 



I/O n-|-2 joiiTi) ^, o, 

3 8m e 



(42) 
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with 



j^(m) = Bmldv v""-' ^\v; m, d*) , 



(43) 



(44) 



and 



jp{m) = Bmjdvv"'+^ ^'^{v;m,d*)E{v;m,d*) . 



(45) 



Except for the factor Bm, which is 



Bm = 



F^,o r(2-f)r(f 



i2 ' 



(46) 



the coefficients i(f,{rn), jaijn)., and jpim) are precisely the integrals denoted 
respectively as Jo,3{'m,d*), J4.3{m,d*), and /i(m,rf*), in I. The quantity Sd = 
2 7r''/^/r((i/2) in Eq. (46) (with d = m, e.g.) means the surface area of a 
(Z-dimensional unit sphere. 

Our two-loop results for the remaining Z-factors, and Z^-, can be written 

as 



n+2 u n+2 

Zj-r-^Jfh = i H 

^ * 3 2e 3 



n-|-5 1 Jui'fn) 
~6~ e2 2^ 



+ 0(^3) (47) 



and 



9 2e 

+ 0(«3), 



/n+8 3 \ 2 5n+22 J„(m) 



V 9 2e 



27 2e 



(48) 



with 



Ju{m) 



^ C^ + ^(2-f) 



(49) 
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where juiiTi) means the integral 



i«M = I dvv-'^\v;m,d*)e{v;m) . (50) 





oo 



For the values m = 2 and m = 6, the above integrals j^, j^r, jp, and ju can be 
computed analytically (cf. I and Appendices B and C). This gives 



4 8 
i0(2) = - , j0(6) = - 



1 - 3 In - 
3 



(51) 



, , 128 , , 448 

M2) = — , j.(6) = , (52) 



Jp(2) = - , Jp(6) = - 



4 

l + 61n- 
3 



and 



(53) 

Jn(2) = -ln^, ^■-(6) = -Q + ln^), (54) 



Ju(2)=ln^, J„(6) = ^-3 1n^. (55) 
3 D 3 



For other values of m we determined these integrals by numerical integration 
in the manner explained in Appendix E. The results are listed in Table 1. 

Table 1 

Numerical values of the integrals j^{m), . . . , Ju{m). 



m 


1 


2 


3 


4 


5 


6 


7 


34> 


1.642(9) 


1.33333 


1.055(6) 


0.803(7) 


0.57(4) 


0.36521 


0.17(4) 




1.339(4) 


4.74074 


10.804(3) 


20.067(7) 


32.95(4) 


49.77778 


70.74(7) 


jp 


0.190(6) 


0.88889 


1.999(9) 


3.464(1) 


5.23(4) 


7.26958 


9.53(6) 


ju 


-0.203(7) 


-0.40547 


-0.624(2) 


-0.880(1) 


-1.21(1) 


-1.72286 


-2.92(4) 




0.383(8) 


0.28768 


0.200(8) 


0.119(8) 


0.04(3) 


-0.02971 


-0.09(9) 



3.3 Beta function and fixed-point value u* 

^From the above results for the renormalization factors the beta function 

(3u{u) = fid^loU (56) 
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and the exponent functions 

Vt = A^^mIo ^^^t, t = 0, cr, p, T, u, 



(57) 



can be calculated in a straightforward manner. Here (9^|g denotes a derivative 
at fixed bare values ctq, tq, po, and uq. Since we employed minimal subtraction 
of poles, the exponent functions satisfy the following simple relationship to 
the residua (Res) of the Z factors: 

r)^{u) = -uduReSe=o[Z^{u)] , c = 4>, a, p, r, u. (58) 



For the beta function, which is related to rju via /9„(m) — —u[e + r)u{u)], we 
obtain 



n+8 3 
~9~ 2 



5n+22 
27 



1 n+2 
^ 24 ~3~ 



8 (m+2) 



J„(m) 



(59) 



Upon solving for the nontrivial zero of we see that the infrared-stable fixed 
point is located at 



9 



,_2e _9_ 8e; 
^ " 3 n+8 ^ 27 Vn+8 
1 n+2 
^ 24 ~Y~ 



[ 5n+22 
3 J^(m) 

jaim) 



8 (m+2) 



+ 0{e') 



(60) 



We refrain from giving the resulting lengthy expressions for the exponent 
functions here. The values r/* = f]^{u*) of these functions at the infrared-stable 
fixed point are presented in Eqs. (61)-(66) below. 



4 Critical exponents 



4-1 Analytic e-expansion results 



The fixed-point value (60) can now be substituted into the exponent functions 
(58) that are implied by our results (40)-(42) and (47) for the renormalization 
factors to obtain the universal quantities rj* — r]i{u*), i = (p, a, p, and r. 
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Recalling how these are related to the critical exponents (cf . I) , one arrives at 
the e expansions 

Vi2 =% = , "i ^ + O e ) , 61) 

Via — 4 — = rrr -, r e + Oie ) , 62 

" 2 + 77* (n + 8)2 2m (m + 2) ^ ^ ' ^ ^ 



1_9- *-_!^ n+2 7n+20 ^ ^ n+2 J^(m)_ 
^ " " n+8 ' 2 (n+8)2 \^ n+8 "^""^ ^ n+8 8 (m+2) 

m (n+2) + 4 (4-n) , J 9 , , 



2 + 77* 1 71+2 e 1 77 + 2 f ^ , 777+20 , , , 

- ' ' ■'77+2 + 4 —Ju{m) 



4 (2 + 7;*) 4 77+8 8 16 (71+8)2 



77+8 



77+2 . 



(m) 



^ 777 77 +2 

4 77+8, 



2777(777+2) 



^e2 + 0(e3) , 
(64) 



and 



V _ * _ ^ + 2 



Ja("0 3jp(rn) ^^(777) 



8777.(777+2) 777 8—777 



e2 + 0(e^).(65) 



The anisotropy exponent 6 — ^1^/^12 of Eq. (2) is given by 



2 + ril 1 77 + 2 



2(77 



+ 



8777 (777 + 2) 8 — 777 



£2 + 0(6^), (66) 



and for the correction-to-scaling exponent, we obtain 



a;,2 = /3>*)=e 



36 e2 577.+22 ^ , , 

3 ^("T.) 



77 + 



27 



1 77+2 

^ 24 ~3~ 



i (m+2) 



- J0(^) 



+ 0(e' 



(67) 



Our rationale for denoting the latter analog of the usual Wegner exponent as 
LJ12 is the following: It governs those corrections to scaling that are weaker 



by a factor of { 



-^12 



|^|f/2t^i2 ^jian ^jig leading infrared singularities. Since 
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^ m2/e .g j^g^^^pg^i introduce also the related correction-to-scaling 
exponent 



<^/4 = 



e 



(68) 



In the case of an isotropic Lifshitz point (cf. Sec. 4.5), in which only the 
correlation length is left, this exponent retains its significance and becomes 
the sole remaining analog of Wegner's exponent. 



4-2 The special cases m — 2 and m — 6 



For these two special cases, a two-loop calculation was performed in Ref. [14]. 
In order to compare its results with ours, we must recall that these authors 
took = Tn — e\\ and d± = 4 — m/2, with m = 2 and m = 6. Accordingly, our e 
must be identified with e = 2 ey. If we substitute the analytic values (51)-(55) 
of the integrals j^, ■ ■ ■ , Ju into our results (61)-(64), the latter reduce to 



77;2(m=2) 



4 n + 2 
9 (n + 8)2 



(69) 



77/4 (m=2) = - 



8 n + 2 2 ^/ 3N 



(70) 



vi2{'m^2) 



1 n+2 

2 ^ 4 (n+8) 
7n+20 , 4 



e + 



2 (n+2) rn^ 



vn+8)3 
+ 0(e=^) , 



131 n 35 
16 ^ 216 ^ 27 



(71) 



216 27 



, 1 n+2 e n+2 115 ri 19 

1/(4 (m=2) = - H h 7 \ ^ 

' 4 n+8 8 (n+8)3 Ll6 

7n+20 , 41 2 ^, 3^ 
+ ^^ln-J e' + 0(e=^). 



(72) 



and 



77/2 (m=6) = 8 Q - In^^ 



(n + 8)^ 



14 n + 2 2 
^M(-=6) = --^^^e^ + 0(e3), 



(73) 
(74) 



16 



ui2{m=6) 



4 (n+8' 



e + 



n+2 ^ ^ 2(n+2) Tn^ 
n+8)3 [16 



+ 



331 n 547 

+ 



144 



72 



23n+88 , 4 

in 

I 3 



(75) 



1 n+2 e n+2 



i'i4{m=6) = - H - + 



4 n+8 8 4(n+8)3 



835 n 1009 



4 108 54 



(19n+56) In- 
3 



+ 0(e^), (76) 



respectively. Our resuits (69), (70), (73), and (74) for tlie correiation exponents 
are consistent witli Eqs. (57), (56), (61), and (60) of Ref. [14]. However, tlie 

terms of tlie correlation- length exponents given in its Eqs. (58), (59), (62), 
and (63) are incompatible with ours. These discrepancies have two causes. 
There is a sign error in Mergulhao and Carneiro's [14] definition (50): The 
term in its second line should be replaced by its negative; only then are their 
general formulas (54) and (55) for the correlation-length exponents uia (= our 
z/;4) and ui/s (= our 1^12) correct. With this correction, these formulas yield 
results in conformity with ours if m = 6. However, in the case m = 2, another 
correction must be made: We believe that in their Eq. (C7) for the integral 
the prefactor of the multiple integral on the right-hand side is too small by a 
factor of 2*^. This entails that the logarithm term ln(16/3) of the pole in 
their Eq. (CIO) gets modified to ln(64/3). With this additional correction, the 

terms following for m = 2 from their corrected Eqs. (54) and (55) turn out 
to be consistent with ours.[3] 



4-3 Numerical values of the second- order expansion coefficients 



In order to analyze further the above results, let us denote the coefficients of 
the terms of the exponents A = z//2, jyu,- . . , (p as C2^\n,m), so that, for 
example, 

= i + ^ ^ 6 + V, m) + 0{^) . (77) 



The numerical values of these coefficients for m = 1, . . . , 6 are listed in Table 
2 for the case n = 1. 



^ We are grateful to C. E. I. Carneiro who checked the calculations of Ref. [14] and 
confirmed the correctness of our results. 
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Table 2 

Numerical values of the second-order expansion coefficients C2^\n=l,m) 



m 





1 


2 


3 


4 


5 


6 




0.043 


0.037113 


0.032158 


0.027744 


0.023677 


0.01987 


0.01626 




— 


0.015867 


0.013335 


0.011083 


0.009010 


0.00707 


0.00524 




-0.090 


-0.062429 


-0.039810 


-0.019277 


-0.000059 


0.01817 


0.03564 




0.006 


-0.00155 


-0.008137 


-0.014196 


-0.019926 


-0.02541 


-0.03070 




0.077 


0.065533 


0.056085 


0.047668 


0.039911 


0.03265 


0.02576 






0.023210 


0.004723 


-0.011535 


-0.026221 


-0.03965 


-0.05203 




-0.630 


-0.482459 


-0.361628 


-0.253233 


-0.152736 


-0.05818 


0.03190 



In deriving the coefficients of the critical exponents and 7^, we utilized 

the familiar hyperscahng and scahng relations [28,7] 

ai = 2- {d-m)ui2-m , (78) 

and 

7i = 1^12 (2 - rii2) = m (4 - riu) , (80) 

respectively. Specifically, the first one, Eq. (78), in conjunction with Eqs. (63) 
and (64) yields 

a, = ^ i + Ct^\n, m) + 0{e') , (81) 
n+Q 2 

with 

C^"'^ ^) ^ ^ ^ ~ ~ ^'"''^ ^2"''^ {n, m) . (82) 

Note that, to first order in e, the expansions of the critical exponents are 
independent of m. This means that one can set m = 0. Hence the expansions 
to first order in e of all critical exponents of the Lifshitz point that have well- 
defined analogs for the usual isotropic (m — 0) critical theory coincide with 
those of the latter, which can be looked up in textbooks [36]. Examples of 
such critical exponents are Vi2i o;;, and 7^. 

The source of this m independence is the following. The operator product ex- 
pansion (OPE) of the theory considered here, for e > 0, is a straightforward 
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extension of the familiar one of the isotropic {m — 0) 0^ theory. Proceeding 
by analogy with chapter 5.5 of Ref. [37]), one can convince oneself that the 
0(e) corrections to the critical exponents are given by simple ratios of OPE 
expansion coefficients. These do not require the explicit computation of Feyn- 
man graphs but follow essentially from combinatorics. For critical exponents 
with an m = analog, this has the above-mentioned consequence. The dif- 
ference between the cases of a Lifshitz point and of a critical point manifests 
itself in the 0(e) expressions of these exponents only through the modified, 
TO-dependent value of e — 4 — d + m/2, a difference that disappears for m — 0. 

In Fig. 1 the coefficients €2^^ of the exponents X — U12, ai, (3i, and 7; for the 
case n = 1 are displayed as functions of m. The results indicate that these 
coefficients depend in a smooth and monotonic fashion on m, approaching the 
familiar isotropic m = values linearly in m. Owing to the above-mentioned 
independence of the 0(e) expressions of these critical exponents, this behavior 
carries over to the numerical estimates one gets for the critical exponents 
in three dimensions by extrapolation of our O(e^) results. We shall see this 
explicitly shortly (see Sec. 4.6). However, before turning to this matter, let 
us briefly convince ourselves that our analytic two-loop series expressions for 
those critical exponents, renormalization factors, etc. that remain meaningful 
for m = go over into the corresponding well-known results of the usual 
isotropic 0^ theory for a critical point. 



.075f 

0.05 
0.025 



-0.025 
-0 . 05 
-0 .075 



▲ 
♦ 



▲ 



★ 
▲ 
♦ 



▲ 
♦ 



3 
I 



m 



Fig. 1. Second-order coefficients C2'^\n,m) of the exponents A = 1/12 (♦), ai {'k), 
Pi (■), and 7; (a) for n = 1 and m = 0, . . . , 7. 
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4-4 The limit m — > 

In this hmit, d 'perpendicular' coordinates, but no 'parallel' ones remain. 
Hence x can be identified with x±, and the free propagator (7) becomes 

G{x) = x-^+' $(z;=0; m=0, d=4 - e) = | t^'^'"" . (83) 

The V = {) value of ^{v) can be read off from the A; = term of the Taylor 
series (12). The expression on the far right of Eq. (83) is indeed the familiar 
massless free propagator [—A]~^{x). 

We must now determine the m — > limits of the integrals j^, j^, jp, and Ju 
(i.e., ju), in terms of which our analytic results given in Sees. 3.2, 3.3, and 4.1 
are expressed. We assert that the correct limiting values are 

j^(0) = 2, j,(0)=j,(0)=i„(0)=0, J„(0) = ^. (84) 

To see this, note that m-dimensional integrals / d™t; fm{v) = Sm Jo^dv fm{v) 
should approach their zero-dimcnsional analog, namely fo{v=0), as m — >■ 0. In 
the case of j^, we have /o(0) = 5*4 $^(0; 0, 4)/Fq2q = 2. The respective values 
/o(0) for ja- and jp vanish because of the exphcit factors of and v"^ appearing 
in their integrands. On the other hand, the vanishing of ju(0) is due to the 
factor Q[v) of its integrand and the fact that 0(0) = according to the Taylor 
series (28). Finally, the value of Ju(0) given above follows from ju(0) = via 
Eq. (49). 

If one sets m = in the series expansions of quantities whose analogs retain 
their significance in the case of the usual critical-point theory, e.g., in Eqs. (40), 
(47), (48), (59), (60), (63), and (67) for Z^, Z^Z^, Z^ZlZf^ u\ r^;, and 
ujn-i utilizing the above values (84) of the integrals, one recovers the familiar 
two-loop results for the standard (f)^ theory. 

Let us also mention that (because of the factor Sm ~ m in 3^) the integrals 
Jo-, ip, and iu vanish linearly in m as m — > 0. Therefore, quantities like Z„Z^^ 
ZpZ^Zy^ ^ Tju, uii, T]*, 9, etc. that involve ratios such as ja{m)/m have a finite 
m — > limit. It is conceivable that the m — > limits of these quantities might 
turn out to have significance for appropriate problems. However, we shall not 
further consider this issue here. 
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4-5 The case of the isotropic Lifshitz point 



In the case of an isotropic Lifshitz point one has d = d\\ and d± = 0. In a 
conventional e expansion one would expand about d*{8) = 8, setting d = d\^ = 
8 — e||. Results to order e| that have been obtained in this fashion for the 
correlation exponents r]i2 and rji^ and the correlation-length exponents 1^12 and 
1/14 can be found in Ref. [19] 

Since the constraint d = d\\ implies that both d and (iy vary as ey is varied, it 
may not be immediately clear that results for this case can be extracted from 
our e expansion at fixed d\^ = m. Let us choose a fixed m = 8 — ey and utilize 
the e expansion. This yields 

A(n,m,d) = A(°) + A(^)He + A(2)(^,^)e' + C'(e3) (85) 



with e = 8 — (i— e||/2, where A means any of the critical exponents considered in 
Sec. 4.1 that remain meaningful in the case of an isotropic Lifshitz point, such 
as 77/4, t'i4, V', «z, Pq-i and cum- As indicated in Eq. (85), the coefficients of the 
terms of orders e° and e do not depend on m. We now set d = m, which implies 
that e = e||/2. The upshot is the following: In order to obtain from our e- 
expansion results the dimensionality expansion of the critical exponents of the 
isotropic Lifshitz point about li = 8 to second order, we must simply replace 
the second-order coefficients X^'^\n,m) by their limiting values z/(2)(n, 8- ") and 
identify e with e||/2 = (8 - d)/2. 

The hmiting values of the integrals j^p, ja, jp, and Ju are 

i^(8-) = 0, j,(8-) = 96, jp(8-) = 12, J„(8-) = -l. (86) 



To see this, note that the factor Bm appearing in j^,. . . , ju varies ~ (8— m) 
near m = 8. In the case of the integral that multiplies Bm has a finite 
m 8 limit, so j<^(8^) vanishes. By contrast, the corresponding integrals 
pertaining to ja and jp have a pole of first order at m = 8. We have in jo-. 



dv 1;™+^ ^%v; m, d*{m)) = I dw w''"^ f{w) [1 + 0(8 - d)] 



7—m 



m 



+0 



m 



57) 



^ We have checked these results by means of an independent calculation, using 
dimensional regularization and minimal subtraction of poles. 
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where f{w) is the function f{w) = [w ^ 8, 8)]^ whose value /(O) = 

(2 7r)-i2 follows from Eq. (22). Using this together with B'{8) = -3 x 2^^ n^^ 
yields the above result for ja(8~). The value of jp(8~) follows in a completely 
analogous manner. 

The computation of J„(8^) is somewhat more involved because the integral 
giving ju/ B„i has poles of second and first order. The second-order pole, due 
to the appearance of the term ~ In v in the large-v form (30) of the scaling 
function 6(1;), produces a first-order pole in j^, which cancels the pole result- 
ing from the contribution — | ^(2—^) to J„ [cf. Eq. (49)]. The first-order pole 
of the integral results in the finite value (86) of Ju(8~). 

Upon substituting the values (86) into the respective e-expansion results of 
Sec. 4.1 and setting 2e — e\\ — S — d,we recover indeed the results of Ref. [19] 
for the correlation exponent 77^4 and the correlation- length exponent 

3 n + 2 

Mm=d) = -^7—^ ej + 0(ej|) , e,, = 8 - d , (88) 



and 



Pu{m=d) 



1 n + 2 
4 ^ 16(n+8) 



(n + 2)(15n2 + 89n + 4) 
960(n+8)3 



(89) 



Furthermore, we can infer the previously unknown series expansions of the 
remaining exponents of the isotropic Lifshitz point. Specifically for the wave- 
vector exponent, we find that 



Another significant exponent is the crossover exponent if. Its 8 — expansion 
follows from Eq. (90) via the scaling law ip = I'u/ Pg. The one of the correction- 
to-scaling exponent (68) becomes 

, 202+41 n 2 ^/ 3x 

uuim^d) = 6|| + el + 0(6}) . (91) 



4.6 Series estimates of the critical exponents for d = 3 dimensions 



We now wish to exploit our e-expansion results of the foregoing subsections 
to obtain numerical values of the critical exponents in d = 3 dimensions. We 
shall mainly consider the cases of uniaxial Lifshitz points (m = 1) for order- 
parameter dimensions n — 1,2,3 and of biaxial (m = 2) Lifshitz points for 
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n = 1. Of particular interest is the case m = n = 1, which is reahzed by 
the Lifshitz point of the ANNNI model [7,38] and is encountered in many 
experimental systems. 

Cases with m > 2 are of limited interest whenever n > 2, for the following 
reason. If a Lifshitz point exists, then low-temperature spin-wave-type exci- 
tations whose frequencies vary as Uq = aoqj + q'j_ as q — > must occur. By 
analogy with the Mermin- Wagner theorem [39] one concludes that such excita- 
tions would destabilize an ordered phase in dimensions d < d^:{m) = 2 + m/2, 
ruling out the possibility of a spontaneous breaking of the 0{n) symmetry 
at temperatures T > for such values of d. (This conclusion is in complete 
accordance with Grest and Sak's work [40] based on nonlinear sigma models.) 
Hence in three dimensions one is left with the case m = 1 of a uniaxial Lifshitz 
point if 72 > 2. 

In Table 3 we list numerical estimates of the critical exponents for d = 3, 
n = 1, and m = 1, 2, . . . , 6. For comparison, we also included the m = values 
of those critical and correction-to-scaling exponents that go over into their 
standard counterparts u, 7, a, P, and u for a critical point.0 As is explained 
in the caption, these estimates were either obtained by setting e = d* (m) — 3 
in the O(e^) expressions of the exponents or else via [1/1] Fade approximants. 

According to Table 2, the coefficients of most of these series with n = 1 
do not alternate in sign. Exceptions are the ones of ai and UJ12 for small 
values of m, that of (3i for m = (i.e., of the usual critical index (3), and the 
one of !f for larger values of m. For d = 3, the second-order contributions 
grow very rapidly as m increases because of the factor = (1 -|- m/2)^. 
Therefore the numerical estimates become less reliable for large m. This effect 
is more pronounced for [1/1] estimates from non-alternating series than for 
the corresponding direct evaluations at = 3 (marked by superscripts (e)). 
The better-behaved expansions yield smaller differences between these two 
kinds of estimates. In unfavorable cases with rather large e we reject the [1/1] 
estimates for the non-alternating series, which tend to overestimate the values 
of the corresponding exponents. Instead we prefer the direct evaluations at 
d = 3. 

A reversed situation occurs for ai and UJ12 with m = 0, 1,2. The respective 
series are alternating; they have negative O(e^) corrections, which tend to 
underestimate the values of the exponents for d = 3 in direct evaluations of 
the O(e^) expressions. On the other hand, the [1/1] approximants for these 
seriesQ seem to do a better job, suppressing the influence of the second-order 

^ Numerical estimates of the correlation exponents r]i2 and rjn are not included in 
Table 3, as they can be found in I. 

^ The e expansions of ai and UJ12 start at order e. We add unity to these series, con- 
struct the [1/1] approximants, and subsequently subtract unity from the resulting 
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Table 3 

Numerical estimates for the critical exponents with n = 1 and d = 3. The values 
marked by superscripts (e) were obtained by setting e = l + m/2in the expansions 
to order of the exponents; those marked by superscripts [1/1] were determined 
from [1/1] Fade approximants whose parameters were fixed by the requirement that 
the respective expansions to second order in e are reproduced. 



m 









1 




2 




3 




4 


5 




6 


(e) 
^12 


0. 


.627 


0, 


.709 


0. 


.795 





.882 





.963 


1.035 


1 


.093 


fl/ll 
'^12 


0. 


.673 


0, 


.877 


1. 


.230 


1 


.742 


2 


.19 


2.26 


2, 


.02 


(e) 
'^14 






0. 


.348 


0. 


.387 





.423 





.456 


0.482 





.500 


n /ii 
lA ' ' 
n4 






0. 


.396 


0. 


.482 





.561 





.606 


0.609 





.585 


(A 

ir 


1. 


.244 


1. 


.397 


1. 


,558 


1 


.715 


1, 


.859 


1.983 


2, 


.08 




1, 


.310 


1, 


.609 


2. 


,02 


2 


.46 


2 


.78 


2.86 


2 


.75 


4^ 


0, 


.077 


0, 


.110 


0. 


,174 





.296 


0, 


.499 


0.806 


1 


.24 




0, 


,108 


0, 


.160 


0. 


,226 


0, 


.323 


0, 


.499 


0.94 


4, 


.6 




0. 


.340 


0, 


.247 


0. 


,134 


-0. 


,005 


-0. 


,18 


-0.39 


-0. 


,7 




0. 


.339 


0. 


.246 


0. 


,131 


-0. 


,029 


-0. 


,28 


-0.75 


-2. 





<^(^) 






0. 


.677 


0. 


,686 





.636 





.514 


0.306 





.001 








0. 


.715 


0. 


,688 





.654 





.628 


0.609 





.595 


^12 


0. 


.370 


0, 


.414 


0. 


,553 





.240 


-0. 


,255 


-0.930 


-1. 


,786 


,,[1/1] 

^12 


0. 


.614 


0, 


.870 


1. 


,161 


1 


.313 


1 


.439 


1.545 


1 


.635 



corrections in a correct way. We believe that ci;| ' \m=n=l) ~ 0.160 belongs 
to our best numerical d — 3 estimates that are obtainable from the individual 
e expansions. 

In the case of a;;2, the second-order correction is much larger. While therefore 
less accurate numerical estimates must be expected, the structure of the e 
expansion for uji2 suggests nevertheless that this correction-to-scaling exponent 
should have a larger value than its m = counterpart uj for the critical point. 
(Recall that the latter has a value close to 0.8 [41]). Our best estimate is the 
[1/1] value 

oui2 ~ 0.870 . (92) 

In order to obtain improved estimates we proceed as follows. We choose the 
"best" d — 3 estimates we can get from the apparently best-behaved e expan- 

numerical values of the approximants. 



24 



sions of certain exponents, express the remaining critical indices in terms of 
the former, and compute their imphed values. Thus we select from Table 3 
the numbers 

uil\d=3;m=l,n=l) -0M8 and af ^^'(3; 1, 1) ~ 0.160 , (93) 

which we complement by our estimate 

77,4(3;,!,!)^ -0.019 (94) 

from I. Substituting these into the second one of the scaling relations (80) for 
■ji and the hyperscaling relation (78) for ai yields 

7;(3;1,1) ~ 1.399 and z/«2(3; 1, 1) ~ 0.746 , (95) 

respectively, from which in turn the values 

77^2(3; 1, 1) ~ 0.124 and A(3; 1, 1) ~ 0.220 (96) 

follow via the scaling relations 77^2 = 2 — ^1/1^12 and Pi = {2 — ai — li)/2. 
Likewise, the choices 

¥^(3; 1, 1) ~ = 0.677 and z/i4(3; 1, 1) ~ = 0.348 (97) 

give for the wave-vector exponent^ 

= vi^l^ (98) 

the estimate 

/5,(3; 1,1) -0.514, (99) 

which is fairly close to the value /3q ~ 0.519 of I. We consider the values 
(92)-(97) and (99) as our best estimates for these 10 critical exponents. 

Table 4 presents an overview of our numerical findings. For convenience, the 
mean-field results are included along with the values the e expansions to first 
and second order take at e = 3/2. The case m = 1 with n = 2 corresponds to 

Note that Eq. (77) of I, which recalls the conventional definition of contains 
a misprint: the variable r should be replaced by q. 
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Table 4 

Critical exponents for m = 1, n = 1, 2, 3, and d = 3. The row marked MF gives the 
mean-field values; rows marked 0(e) and O(e^) list the values obtained by setting 
e = 3/2 in the expansions to first and second order in e, respectively. The remaining 

rows contain the estimates from [1/1] Fade approximants and our "best" estimates 
('Seal.') obtained via scaling relations in the manner explained in the main text. 









71 


= 1 




71 


= 2 


ft 


= 3 




MF 


0(e) 




fl/ll 


Seal. 


0(e) 










1 

2 


f) fi2B 


n 70Q 


877 


74fi 


65 


757 


67 


7Q8 


IJl A 


1 

4 


f) 31 3 


n 348 


3Qfi 


348 


325 


372 


335 


3Q2 







0.25 


0.110 


0.160 


0.160 


0.15 


-0.047 


0.068 


-0.178 


A 


1 

2 


0.25 


0.247 


0.246 


0.220 


0.275 


0.276 


0.295 


0.301 




1 


1.25 


1.397 


1.609 


1.399 


1.3 


1.495 


1.34 


1.576 




1 

2 


0.625 


0.677 


0.715 


0.677 


0.65 


0.725 


0.67 


0.765 




1 
2 


0.5 


0.519 




0.514 


0.5 


0.521 


0.5 


0.521 


ril2 








0.039 




0.124 





0.042 





0.044 


Vl4 








-0.019 




-0.019 





-0.020 





-0.021 


^72 





1.5 


0.414 


0.870 




1.5 


0.4G6 


1.5 


0.517 



the Lifsliitz point of the axial next-nearest-neighbor XY (ANNNXY) model, 
which Selke studied many years ago by means of Monte Carlo simulations [42] . 

In Table 5 we have gathered the available experimental results for critical ex- 
ponents together with estimates obtained from Monte Carlo calculations and 
high-temperature series analyses. As one sees, our field-theory estimates are 
in a good agreement with the Monte Carlo results. The experimental value 
for ai{2>] 1, 1) deviates appreciably both from all theoretical estimates (includ- 
ing ours) as well as from the Monte Carlo results, and is probably not very 
accurate. On the other hand, the very good agreement of our field-theory es- 
timates with the most recent Monte Carlo estimates by Pleimling and Henkel 
[27] (which we expect to be the most accurate ones) is quite encouraging. Cer- 
tainly, renewed experimental efforts for determining the values of the critical 
exponents in a more complete and more precise way would be most welcome. 



5 Concluding remcirks 

The field-theory models (3) were introduced more than 25 years ago to de- 
scribe the universal critical behavior at m-axial Lifshitz points [19]. While 
some field-theoretic studies based on the e expansion about the upper critical 
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Table 5 

Values of critical exponents for uniaxial Lifshitz points (m=l). 'Exp' means summa- 
rized experimental results, taken from Ref. [43]; 'HT' denotes the high-temperature 
series estimates of Ref. [44]. The rest are the Monte Carlo results of Refs. [45] (MCI), 
[46] (MC2), [27] (MC3), and [42] (MC), respectively. 







ai 


A 


11 










n = 


1 




Exp 




0.4-0.5 






0.60-0.64 0.44-0.49 


HT 


0.41±0.03 


0.20±0.15 




1.62±0.12 


0.5 


MCI 






0.21±0.03 


1.36±0.005 




MC2 


0.33±0.03 


0.2 


0.19±0.02 


1.4±0.06 




MC3 




0.18±0.03 


0.235±0.005 


1.36±0.03 










n = 


2 




MC 




O.libO.14 


0.20±0.02 


1.5±0.1 





dimension d*{m) emerged soon afterwards, these were limited to first order in 
e, or restricted to special values of m or to a subset of critical exponents, or 
challenged by discrepant results (see the references cited in the introduction) . 
Two-loop calculations for general values of m appeared to be hardly feasible 
because of the severe calculational difficulties that must be overcome. 



Complementing our previous work in I, we have presented here a full two- 
loop RG calculation for the models (3) in d — d*{m) — e dimensions, for 
general values of m e (0, 8). This enabled us to compute the e expansions of 
all critical indices of the considered m-axial Lifshitz points to second order 
in e. We employed these results in turn to determine field-theory estimates 
for the values of these critical exponents in three dimensions. Although the 
accuracy of these estimates clearly is not competitive with the impressive 
precision that has been achieved by the best field-theory estimates for critical 
exponents of conventional critical points (based on perturbation expansions 
to much higher orders and powerful resummation techniques [41,4]), they are 
in very good agreement with recent Monte Carlo results for the uniaxial scalar 
case m — n — 1. We hope that our present work will stimulate new efforts, 
both by experimentalists and theorists, to investigate the critical behavior at 
Lifshitz points. 

There is a number of promising directions in which our work could be ex- 
tended. For example, building on it, one could compute other universal quan- 
tities, such as amplitude ratios and scaling functions, via the e expansion. 

A particular interesting and challenging question is whether the generalized in- 
variance found by Henkel [25] for systems whose anisotropy exponents take the 
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rational values ^ = 2/p, p = 1, 2, 3 . . . can be generalized to other, irrational 
values. That such an extension exists, is not at all clear since the condition 
6* = 2/p is utilized in Henkel's work to ensure that the algebra closes. But if 
such an extension can be found, then the invariance under this larger group 
of transformations should manifest itself through properties of the theories' 
scaling functions in d < d*{m) dimensions, which could be checked by means 
of the e expansion. Furthermore, even if an extension cannot be found, one 
should be able to benefit from the invariance properties of the free theory 
(with 6 — 1/2) when computing the e expansion of anomalous dimensions 
of composite operators in a similar extensive fashion as in the case of the 
standard critical-point 0^ theory [47-49] . 

An important issue awaiting clarification arises when m > 2. In the class of 
models (3) studied here, the quadratic fourth-order derivative term was taken 
to be isotropic in the subspace M™. However, in general further fourth-order 
derivatives cannot be excluded. That is, the term {(Tq/2) (A||0)^ should be 
generalized to 

((Jo/2) Wa {didj(j>)dkdi (j> , (100) 

where the summation over a comprises all totally symmetric fourth-rank ten- 
sors Tj^^li compatible with the symmetry of the considered microscopic (or 
mesoscopic) model. The isotropic fourth-order derivative term corresponds to 
Ti^ki^^ = (% 6ki + 6ik Sji + 6ii 6jk)/3 with wi = 1. 

In order to give a simple example of a system involving a further quadratic 
fourth-order derivative term, let us consider an m-axial modification of the 
familiar uniaxial ANNNI model [38,7] that has competing nearest-neighbor 
(nn) and next-nearest-neighbor interactions along m equivalent of the d hy- 
percubic axes (and only the usual nn bonds along the remaining d — m ones). 
Owing to the Hamiltonian's hypercubic (rather than isotropic) symmetry in 
the m-dimensional subspace, just one other fourth-order derivative term, cor- 
responding to the tensor Tj^jli = Sij5kiSii, must generically occur besides the 
isotropic one, in a coarse-grained description. The associated interaction con- 
stant W2 is dimensionless and hence marginal at the Gaussian fixed point. 
To find out whether the nontrivial {u* > 0, W2 — 0) fixed point considered 
throughout this work remains infrared stable, one must compute the anoma- 
lous dimension of the additional scaling operator that can be formed from 
the above two fourth-order derivative terms. This issue will be taken up in a 
forthcoming joint paper with R. K. P. Zia [50], where we shall show that the 
associated crossover exponent, to order e^, is indeed positive. Hence, devia- 
tions from W2 = correspond to a relevant perturbation at the W2 — 0, u* > 
fixed point, which should destabilize it unless m — 1. 

Finally, let us mention that the potential of the e-expansion results presented 
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in this paper certainly has not fully been exploited here. When estimating 
the values of the critical exponents for = 3 dimensions, we utilized only 
their e expansions for a fixed integer number of the parameter m. However, 
our results hold also for noninteger values of m. Making use of this fact, one 
should be able to extrapolate to points {d — 3,m — integer) of interest in a 
more flexible fashion, starting from any point on the critical curve d = d*{m) 
and going along directions not perpendicular to the m axis. By exploiting this 
flexibility one should be able to improve the accuracy of the estimates. 

Last but not least, let us briefly mention where the interesting reader can find 
information about experimental results. Earlier experimental work is discussed 
in Hornreich's and Selke's review articles [6,7]. A more recent summary of 
experimental results for the critical exponents and other universal quantities 
of the Lifshitz point in MnP and Mno.9Coo.1P can be found in Ref. [43] and 
its references. (These results were partly quoted in Table 5.) However, the 
variety of experimental systems having (or believed to have) Lifshitz points 
is very rich, ranging from ferromagnetic and ferroelectric systems to polymer 
mixtures. A complete survey of the published experimental results on Lifshitz 
points is beyond the scope of the present article. 
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A Series representation and asymptotic expansion of ^{v) 
^Prom the integral representation (8) of the scaling function ^{v) we find 

qy <1± II Q\\ 9± II 

where ey = v/v is an arbitrary unit m- vector. The second form follows via 
rescaling of the momenta; it lends itself for studying the large-v behavior. The 
first one is appropriate for deriving the small- v expansion. 
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Upon utilizing the Schwinger representation 



1 



oo 

y'dse-^^'^W) (A.2) 



for the momentum-space propagator in Eq. (A.l), we can perform the inte- 
gration over to obtain 



oo 



<^{v) = (4 7r)-'^^/2 v-^+^' J ds s-"^!'' j e-^«r(^^"')"' e^n-w . (A.3) 



Doing the angular integrations gives 

oo 

j g-^^e^'^ii-'^li = (2 7:)-'^/^ Jdqq""/^ Jru^iq) e"^"' . (A.4) 



We insert this into Eq. (A.3), expand the exponential c i/C^*'^'') jj^ powers of 
v~'^, and integrate the resulting series termwise over s. This yields 



.,)..-=-™.-^(|)^-|(^..(.)- 



with 







2r(2 + 2A;-e) 



(A.6) 



which is the asymptotic expansion (15). The Taylor series (12) can be derived 
along similar lines, starting from the first form of the integral representation 
(A.l). 



B Laurent expansion of required vertex functions 



In this appendix we gather our results on the Laurent expansions of those ver- 
tex functions whose pole terms determine the required renormalization factors. 
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It is understood that tq and po are set to their critical values Tc = pc = 0. For 
notational simplicity, we introduce the dimensionless bare coupling constant 

uo = a^"'^'^ uo = ZuU (B.l) 



and specialize to the scalar case n = 1. The generalization to the n-component 
case involves the usual tensorial factors and contractions of the standard \(p\^ 
theory and should be obvious. 

We use the notation r{{qj}) for the Fourier transforms of vertex functions 
r{{xj}) (with the momentum-conserving S function taken out): 

r({x,})= / f({g^.})(27r)^5(Eq^.)e^S.-i'^^-^ . (B.2) 

B.l Two-point vertex functions F^^^ and T^^y^^^y 

/^From our results obtained in I we finely 

r^'\q) = ^oqt + ql- + 0(4) (B.3) 



with 



6e 



lQm(m+2) 2 (8 - m) 



+ 0(e° 



(B.4) 



and 

f Q = 0) = ql- + 0{ul) (B.5) 



with 



where Q = is the momentum of the inserted operator (V||0)^/2; i.e., the 
insertion considered is / d'^x(V|| 0)^/2. 

^ We suppress diagrams involving the one-loop (sub)graph since the latter 
vanishes for r = if dimensional regularization is employed, as we do throughout 
this paper. 
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B.2 Four-point vertex function T^'^^ 



The four-point vertex function was computed only to one-loop order in I. To 
the order of two loops it reads ^ 



= «o - ( >co< + 2 p) - ( >c:xo< + 2 P) 



(ii)=(12),(23),(24) L 



Or/- \ "0 r2/- 

Y/2(q.,)-^/2(9^,. 

+ f [/4(qi2,g3)+5P]}. 



(B.7) 



(B.8) 



Here V means permutations (of the external legs). The hatted momenta are 
dimensionless ones defined via 



- \ / 1/4 -1/2 -1 



(B.9) 



and q^j = + qj. The integrals I2 and J4 are given by 



l2iQ) = / 7 rj^ 7 (B.IO) 



and 

^ g± + ?ll {q + Q)^ + iq + Q)\\ 



The pole term of hiQ) can be read off from Eqs. (24) and (89) of I. However, 
in our two-loop calculation I2 also occurs as a divergent subintegral. To check 
that the associated pole terms are canceled by contributions involving one- 
loop counterterms, we also need the finite part of 12- The calculation simplifies 
considerably if the momentum Q is chosen to have a perpendicular component 
only, so that Q = Q e± (which is sufficient for our purposes) .[^ For such values 

^ Previously e± denoted a fixed arbitrary unit d — m vector. For convenience, we 
use here and below the same symbol for the associated d vector whose projection 
onto the perpendicular subspace yields the former while its m parallel components 
vanish. 
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of Q, the integral l2{Q) can be analytically calculated in a straightforward 
fashion, either by going back to Eq. (16) and (19) of I and computing the 
Fourier transform of these distributions, or directly in momentum space, as 
we prefer to do here. For dimensional reasons, we have 



l2{Qe^)^Q-'h{e^) . 



(B.12) 



The integral on the right-hand side is precisely the one written as -F^,e/e in 
Eq. (38). Utilizing a familiar method due to Feynman for folding two denom- 
inators into one (Eq. (A8-1) of Ref. [36]), one is led to 



q(\ + ql + s qj_ ■ e± + s e]_ 



q 



= 2-^+^+' 7r-2+T+2 r ■ / ds 



(47r) 



ql + s{l - s) 



m+2 e 
4 



911 



-./2r(f)r(i)r^(i-i) 



2r(2-e)r(f) ' 

from which the result (39) for g follows at once. 



(B.13) 



The calculation of h{Q, K) is more involved; it is described in Appendix C, 
giving 



+ J„(m) +0(e) 



(B.14) 



where Ju{Tn) is the quantity defined in Eq. (49). 



B. 3 Vertex function F 



(2) 

,2 



Next, we turn to the vertex function r^2 with an insertion of |(0^)g = 
|/ d'^x0^(£c) e^Q'^. To two- loop order it is given by^ 



rl?(9;Q) = i- Q - S - +0{ul) 



+ 0{ul), (B.15) 
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where the hatted momenta are again dimensionless ones, defined by analogy 
with Eq. (B.9). 



C Laurent expansion of the two-loop integral I4 



As can be seen from Eq. (B.ll), the integral l4:{Qi2yQ3) associated with the 
graph involves the divergent subintegral l2{q — qs)- The latter has a 
momentum-independent pole term ~ [cf. Eqs. (B.13) and (39)]. Further- 
more, the graph that results upon contraction of this subgraph to a point 
[which itself is proportional to -^2 (^12)] has contributions of order e° that de- 
pend on qi2. Taken together, these observations tell us that the pole term 
oc of h{qi2,Q3) depends on qu but not on q^^^ Since the pole of 
h{Qi2: Q3) is momentum independent, we can set qs = when calculating the 
pole part of this integral. 

To further simplify the calculation, we can choose to have vanishing parallel 
component again, setting q-^2 = Q The integral to be calculated thus 
becomes 

h{Qer,0) = Jd''xJd''yG{y)G{x-y) e^'^^^-y G'{x), (C.l) 



where G{y) now means the free propagator (5) with ctq = 1. Let us substitute 
the free propagators of the factor G'^{x) by their scaling form (7) and rewrite 
the Fourier integral / d'^y ... as a momentum-space integral, employing the 
Schwinger representation (A. 2) for both of the two free propagators in mo- 

— 1/2 

mentum space. Making the change of variables x\\ ^ v = x^\ Xj_ , we arrive 
at 



00 00 



h{Qe^,0) = Jd'''v<l>\v)Jd'^-'^x^xl"^^ V^^/^^ 



X f Q-qiis+t)+q_L-iix^~2tQe^)~tQ'^ f-q^^ {s+t)+i q^yu^/xT 



9j_ 



Now the momentum integrations /^^ and J^^^ are decoupled and can be per- 
formed in a straightforward fashion. That the latter integral takes such a 
simple form is due to our choice of Q with Qy = 0. Performing the angular 
integrations yields 



It is precisely this gi2-dependent pole term that gets canceled by subtracting from 
the divergent subgraph its pole part. 
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(27r)-^ J dq\\ 



T -<1« / 2 



x±) 4 J „x-2 V xY'^^ . (C.3) 



We replace the Bessel function in Eq. (C.3) by its familiar Taylor expansion 



(C.4) 



integrate term by term over , employing 



(C.5) 



and simplify the resulting ratio of F-functions by means of the well-known 
duphcation formula (6.1.18) of Ref. [35]. This gives 



/ 



m+2 k 
2 



. (C.6) 



The integration over q^^ in Eq. (C.2) is Gaussian. Upon substituting the result 
together with the above equations into (C.2), we get 



oo oo 



h{Qe±,0) = 2'-^-"'7r^^ ld"'v^\v) j ds j dt {s + t)— e 



,2 



(C.7) 



with 



C,(g;s,t)= d^-^'x^x 



exp 



a;^ + 4 i t Q e_L • cc_L 
4(s + i) 



(C.^ 



We first perform the angular integrations and subsequently the radial integra- 
tion of the latter integral, obtaining 



.S + ty 



35 



where we have introduced 

d — m ^ m e . „s 

= = 2 - - - - . (C.IO) 

Next we insert this result into expression (C.7) for I4, and make a change 
of variable s ^ z — s/t The ^-integration then becomes straightforward (see 
Eq. (2.22.3.1) of Ref. [51]), and we find that 



h{Q e±,0) =Q 



re 



-m— 2 e 



/oo 
(Tv ^\v; m, d) ^ ^(m, e) {-v'^f , (C.ll) 
fe=0 

with 



(C.12) 

Owing to the overall factor r(e) and the additional factor r(e/2) of the coeffi- 
cient A^im^ e), the k — ^ term of the above series contributes poles of second 

and first order in e to /4. The remaining terms with A; > 1 yield poles of first 
order in e. Consider first the A; = term. The value of the integral over v may 
be gleaned from I [cf. its Eqs. (3.16) and (4.47)]: 



/ dTv $2(i;; m, d) = ^^-^ ^ — ■ (C.13) 

J ^ ' ' ^ r(f)r(2-e) ^ ^ 



The integral over z in Ao{m, e) can be evaluated explicitly by means of Math- 
ematica [52]. Alternatively, one can change to the integration variable ( — 1/z 
and look up the transformed integral in Eq. (2.21.1.15) of the integral tables 
[51]. The result has a simple expansion to order e, giving 



Ao{m,e) 



r(6/2) 
r(^) 



l + 2e + 0(e2) 



(C.14) 



upon substitution into Eq. (C.12). 



36 



Turning to the contributions with /c > 1, we note that both the scahng function 
^[v;m,d) and the coefficients Ak{m,e) may be taken at d = d* (i.e. e = 0). 
Then the integral J^dz . . . reduces to one and the series J2k'=i becomes the 
function Q{v;m) introduced in Eq. (28). It follows that 

/ d"^v <^'iv; m, d) f: A,(m, e) {-v'f = ^ [1 + 0{e)] , (C.15) 



where ju(m) and Bm are the integral (50) and the coefficient (46), respectively. 
Combining the above results and expanding the prefactors of the integral in 
Eq. (C.ll), we finally obtain the result stated in Eq. (B.14). 



D Asymptotic behavior of Q{v) 



Upon differentiating the series (28) of Q{v;m) termwise and comparing with 
the Taylor expansion (12) of the scaling function $, one sees that the following 
relation holds: 



dQ{v;m) 4 



dv 



V 



^v;m,d*) -^0;m,d*) 



(D.l) 



^From Eq. (12) we can read off the value ^{0;m,d*) = l/r[(m + 2)/4]. Let us 
substitute the asymptotic expansion Eq. (15) of $(t;; m, rf*) into this equation 
and integrate. This yields 

e(vm) ^ -4 1nt^ + ge(m) - 2(-l)^r(2fc) fvy^^ 



The terms of orders v ^ and v ^ agree with those of the asymptotic form (30) 
of 0(i;; m). Hence it remains to show that the integration constant Cq is given 
by 

f rr, _|_ 2\ 

Ce(m) = V \-^) - C^; + In 16 . (D.3) 



To this end an integral representation of Q{v\m) is helpful. Consider the 
integral 



Jq{v- m, e) = 2^+™ t^"^^ \ ^ = V ^ ^ ^ , (D.4) 
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in terms of which Q{v;m.) can be written as 



Q(v: m) = hm 



Je{v;m,€) 



£(ztl 

r(^) 



(D.5) 



and whose large-v form 



2^-^^r(-e) 

r(^) 



V 



2e 



1 + 0(0 



(D.6) 



is easily derived. Insertion of the latter result into Eq. (D.5) gives the value 
(D.3) of Ce. 



E Numerical integration 



The quantities j(j){m), jaijn)-, 3pijn)-i and juiin) in terms of which we expressed 
the series expansion coefficients of the renormalization factors and the critical 
exponents arc integrals of the form J^dvf{v) [cf. Eqs. (43)-(45) and (50)]. 
Their integrands, /, while integrable and decaying to zero asv oo, in general 
involve differences of generalized hypergeometric functions, i.e., differences of 
functions that grow exponentially as f — > oo. Therefore standard numerical 
integration procedures run into problems when the upper integration limit 
becomes large. 

To overcome this difficulty, we proceed in a similar manner as in I. ^Prom 
our knowledge of the asymptotic expansions of the functions ^{v;m,d*), 
S(-u; m, d*), and Q{v;m) wc can determine that of the integrand. Let f^\v) 
be the asymptotic expansion of f{v) to order . Then we have 



oo 



(E.l) 



k=M+l 



We split the integrand as 

oo VQ VQ 

jdvf{v) = jdvfiv)- Jdvf^\v) + Rf\vo), (E.2) 
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where 



(E.3) 



vo 



Then we choose Vq as large as possible, but small enough so that Mathematica 
[52] is still able to evaluate the integral jQ°f{v) dv by numerical integration, 
determine the second term on the right-hand side of Eq. (E.2) by analytical 
integration, and neglect the third one. The asymptotic expansion of the latter 
is easily deduced from Eq. (E.l). It reads 

oo 

<'W„-^E^«o-'. (E.4) 

k=M 



Since the expansion (E.l) is only asymptotic, the value of M must not be 
chosen too large. In practice, we utilized the asymptotic expansions of $, S, 
and © up to the orders v"^^, v~^'^, and explicitly shown in the respective 
Eqs. (22), (25), and (30), and then truncated the resulting expression of the 
integrand f{v) consistently at the largest possible order. As upper integration 
limit Vo of the numerical integration we chose values between 9 and 10. 

As a consequence of the fact that all integrands f{v) have an exphcit factor 
of v"^, the precision of our results decreases as m increases. Furthermore, the 
accuracy is greatest for j^(m), whose integrand's asymptotic expansion starts 
with particularly high power of v ^. The precision is lower for 

ja{fn) and ju{m) because their integrands involve either four more powers of 
V than that of or else the function Q{v) as a factor, whose asymptotic 
expansion starts with a term ~ Inv. 

As a test of our procedure we can compare the numerical values of the integrals 
it produces for m = 2 and m = 6 with the analytically known exact results 
(51)-(54). The agreement one finds is very impressive: Nine decimal digits 
of the exact results are reproduced (even for m = 6) when the numerical 
integration is done by means of the Mathematica [52] routine 'Nintegrate' with 
the option 'WorkingPrecision=40'. However, we must not forget that the cases 
m = 2 and m = 6 are special in that the asymptotic expansions of the 
functions S, and — and hence those of the integrands — vanish or truncate 
after the first term. Hence it would be too optimistic to expect such extremely 
accurate results for other values of m. In the worst cases (e.g., that of jo- (7) 
and ju(7)), the fourth decimal digit typically changes if v is varied in the 
range 9 ... 11. Therefore we are confident that the first two decimal digits of 
the m — 1 values given in Table 1 are correct. For smaller values of m the 
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precision is greater. 
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